Teleportation fidelity as a probe of sub-Planck phase-space structure 



A. J. Scott^B and Carlton M. Caves 2 - 3 -Q 

1 Centre for Quantum Computer Technology and Centre for Quantum Dynamics, 
Griffith University, Brisbane, Queensland 4111, Australia 
2 Department of Physics and Astronomy, MSC07-4220, 
University of New Mexico, Albuquerque, NM 87131-0001, USA 
Department of Physics, University of Queensland, Brisbane, Queensland 4 072, Australia 

We investigate the connection between sub-Planck structure in the Wigner function and the 
output fidelity of continuous-variable teleportation protocols. When the teleporting parties share 
a two-mode squeezed state as an entangled resource, high fidelity in the output state requires a 
squeezing large enough that the smallest sub-Planck structures in an input pure state are teleported 
faithfully. We formulate this relationship, which leads to an explicit relation between the fine-scale 
structure in the Wigner function and large-scale extent of the Wigner function, and we treat specific 
examples, including coherent, number, and random states and states produced by chaotic dynamics. 
| We generalize the pure-state results to teleportation of mixed states. 
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I. INTRODUCTION 
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Quantum teleportation [l|, 0, H, 3 transfers a quantum state from a quantum system owned by Alice to a distant 
system owned by Bob at a cost, per qubit teleported, of one bit of entanglement shared between Alice and Bob 
and two bits of classical communication. Bennett et al. [1] showed that a maximally entangled state shared between 
{ Alice and Bob allows this process, which might naively be presumed impossible under quantum mechanical law: an 
arbitrary quantum state of a system possessed by Alice can be perfectly transferred to a system possessed by Bob 
through local operations and classical communication alone. This apparent violation of the no-cloning theorem 0, Q 
is explained by the complete corruption of Alice's system after teleportation. 

Teleportation has become a fundamental primitive for quantum information processing 0- Although the telepor- 
tation protocol devised by Bennett et al. applies only to finite dimensional systems (e.g., a spin-i particle), it was 
later extended to include continuous- variable systems (e.g., the modes of an electromagnetic field) by Vaidman 2] 
and then by Braunstein and Kimble [3[. In the latter case, a highly squeezed two- mode vacuum state is chosen for the 
entangled resource, and only in the limit of infinite squeezing is the teleported state a perfect replica of the original. 
An appropriate measure of the "quality" of the teleported state is its probability overlap with the original, which is 
called the fidelity. 

Experimental demonstrations [1, [l(J [HI, G2, EH, EH of continuous-variable teleportation have now achieved 
fidelities up to 0.76 ± 0.02 [TJ] when the teleported state is Gaussian. It has been argued that teleportation of 
non-Gaussian states with fidelities above 2/3 is necessary before these experiments no longer afford interpretations 
formulated purely in terms of classical correlations and thus require a quantum-mechanical explanation [l5j |. In 
the present article we show that achieving such high fidelities in teleporting non-Gaussian states requires faithful 
reproduction of the smallest sub-Planck structures in the teleported state's Wigner function. 

The existence and importance of sub-Planck structure in the Wigner function was pointed out by Zurek [p6| . He 
showed that environmental perturbations to a quantum system at the sub-Planck action scale, a ~ 1/^4, where A is the 
areal extent of the phase-space region over which the system state has nonnegligible Wigner function, are enough to 
cause orthogonality between perturbed and unperturbed states and, hence, to drive decoherence. This same sensitivity 
to perturbation has been suggested as a way to improve sensitivity in weak- force detection [It], [H, Gil i m muc h the 
same way that the sub-Planck variation in a squeezed state along one phase-space dimension can be used for this 
purpose. Structures analogous to the sub-Planck structures in the Wigner function of non-Gaussian states have been 
identified in the classical time-frequency domain of an electromagnetic field mode [201 ] , and non-Gaussian states of an 
optical field mode with sub-Planck structure have been generated and investigated by Ourjoumtsev et al. [2ll l22l| . We 
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complement all this previous work by demonstrating the importance of sub-Planck structure to achieving high-fidelity 
continuous- variable teleportation. 

We show that when the teleporting parties share a two-mode squeezed state as an entangled resource, high fidelities 
in teleporting a pure state p require squeezing large enough that the smallest sub-Planck structures in the Wigner 
function of p arc tclcported faithfully. While this connection is reasonable on its face, we make it mathematically 
explicit by showing that the rate of decrease in fidelity is directly related to a natural measure of the fine-scale 
structure in the Wigner function and, reciprocally, to a sensible measure of the large-scale extent of the Wigner 
function, i.e., the size of the region that encompasses all of the nonnegligible support of the Wigner function. This 
explicit connection takes the form 



J dxdp \VWp(x,p)\ 
dxdpW p 2 (x,p) 



\ [(Ax) 2 + (Ap) 2 ] = -J— — = - J / dxdp\\7 W ;(x,p)\ 2 . (1.1) 



dF p 
dt 

In this equation F p (t) is the average fidelity of the teleported state for pure input state p, as a function of a squeezing 
parameter t = 2e~ 2r , where r is the standard squeezing parameter (i.e., t is the twice the ratio of the uncertainty 
in a squeezed quadrature component to the uncertainty in an unsqueezed quadrature); (d F p /dt) t= o is the rate of 
decrease of fidelity away from the perfect fidelity achievable for infinite squeezing (t = 0); (Ax) 2 = (x 2 ) — (x) 2 and 
(Ap) 2 = (p 2 ) — (p) 2 are the variances of x and p ; and W' p (x,p) is the Wigner function of p, normalized to unity with 

respect to the phase-space measure dxdp. The first equality relates (d F p /dt) t= o to a sensible measure of the extent of 
the Wigner function in phase space, and the second equality relates it to a natural measure of the fine-scale structure 
of the Wigner function. The third equality follows from the fact that the obvious measure of the (inverse) area of 
support of a Wigner function has the same value for all pure states, 

J dxd V W l2 (x,p) = i- , (1.2) 

corresponding to one Planck unit of action. That the area of support of the Wigner function is one Planck unit for 
all pure states already tells one that when the large-scale extent is much bigger than a Planck unit, the support must 
be very patchy within that extent, leading to fine-scale structure. 

The explicit connection among teleportation fidelity for pure states, the fine-scale structure of the corresponding 
Wigner function, and the large-scale extent of the Wigner function, all expressed in Eq. (|1.1|) . is the key result of 
this paper. This connection can be generalized to mixed states by using entanglement fidelity [23| as the measure 
of success. With entanglement fidelity in place of fidelity, the first equality in Eq. (|1.1[) still holds, but the second 
and third equalities do not. The third equality can be replaced by a more general measure, which expresses not the 
fine-scale structure of the mixed state being teleported, but rather the fine-scale structure in any purification of that 
mixed state. Indeed, the justification for the generalization to mixed states is to let teleportation fidelity identify for 
us an appropriate measure of the fine-scale structure underlying a mixed state. 

The article is organized as follows. Section |TT] reviews the definitions of ordered characteristic functions and the cor- 
resonding quasiprobability distributions. The general procedure for continuous-variable teleportation and its analysis 
in terms of Wigner functions is outlined in Sec. IIII Al The general analysis is then specialized to squeezed-state tele- 
portation in Sec. IIII Bi and the high-fidelity limit and the requirements on the classical communication are discussed 
in Sec. IIII Cl In Sec. HVI we explore the connection between teleportation fidelity and sub-Planck structure. In Sec. fVl 
we consider specific examples, including teleportation of coherent, number, and random states. Sec. I VII generalizes 
our results to teleportation of mixed states, and we conclude with a brief summary in Sec. IVHl 



II. CHARACTERISTIC FUNCTIONS AND QUASIPROBABILITY DISTRIBUTIONS 

Before proceeding to our analysis of teleportation, we summarize some standard definitions used throughout the 
article — and our particular conventions for those definitions [24| . Consider a mode with annihilation operator a = 
(x+ip) / y/2, where x andp can be thought of as position and momentum operators and are sometimes called quadrature 
components. With this choice, we adopt natural oscillator units in which position and momentum have the same 
(dimcnsionless) units, scaled so that the vacuum uncertainties in position and momentum are Ax = Ap= 1/V2. 

We use the displacement operator 



D(a, a) = exp(ae>J — a* a) , 



(2.1) 
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written as a function of two variables, the first an operator and the second a complex amplitude a = (ai + ict) / y/2. 
The real quantities ol\ and a 2 can be thought of as dimensionlcss c-numbcr variables for position and momentum, 
respectively. We use this decomposition for all Greek variables. The displacement operator generates a coherent state 
with complex amplitude a from the vacuum state, i.e., D(a, a)\0) — \a). 

When the operator slot in the displacement operator is filled by a c-number, we get the two-dimensional Fourier 
expansion function, 

D(/3,a) = e" 13 '- ' 13 = e *( a ^~ a iM , (2.2) 

which satisfies 

d 2 a 



TV 



-D(P,a) =tt6(/3) . (2.3) 



The s-ordered characteristic function for the state p is defined as 

= e W/hr[pD(a,!3)] . (2.4) 

Fors = 0, $ p (a) = ^ p °\a) is the symmetrically ordered characteristic function. The (real-valued) s-ordered quasiprob- 
ability distribution is 1/tt times the Fourier transform of the corresponding characteristic function, 

W«(a) = J (/?)£>(/?, a) = hr[ P D^{a, a)} . (2.5) 

Here the integration measure is d 2 /3 = c?(Re (3) d(lm (3) = dj3\ d/?2/2, and 

D^(a,a) = [ —D^(a,(3)D(f],a) (2.6) 



is the (Hermitian) Fourier transform of the s-ordered displacement operator D^(a,j3) = e s ^^^ 2 D(a,/3). These 
functions satisfy 



d 2 aW { p s) {a) = *W(0) = tip = 1 , (2.7) 
tt/ d'aW^i^W^ia) = j ^$(-)*(^)$W tf) = tr( Pl p 2 ) . (2.8) 



When s = — 1 and s = 0, we have 

D(- 1 )(a,a) = |a)(a|, (2.9) 
£)W(a,a) = 2D(a,a)(-l) a>a D^a,a) = 2 J ^-\a + (3) (a - 0\D{/3,a) . (2.10) 

These allow us to write expressions for the s — and s = — 1 quasidistributions in terms of the coherent-state matrix 
elements of the density operator. These two quasidistributions are called, respectively, the Wigner function and the 
Husimi (or Q) function: 

W p (a) = W^\a) = A j d 2 (3(a-(3\p\a + p)D(P,a) , (2.11) 

Q p (a) = W^Xa) = ~(a\p\a) . (2.12) 

These expressions imply that —2/tt < W p {a) < 2/tt and < Q p (a) < 1/tt. Since the Husimi function is nonnegative, 
it is a probability distribution, rather than just a quasidistribution, albeit one that cannot be too highly peaked. The 
ordered characteristic functions and corresponding quasidistributions were first explored systematically by Cahill and 
Glauber 

When the Wigner function is written as a function of x = a% and p — a?, it is conventional to rescale it by 
a factor of two, i.e., W'{a.\,a.2) = W p (a)/2, so that it is normalized to unity with respect to da±da2, instead of 
d 2 a = da>i da 2 /2. This rescaled Wigner function appears in Eqs. (jl.ip and (|1.2|) and the surrounding discussion, but 
we make no further use of it in this paper. 
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III. CONTINUOUS- VARIABLE TELEPORTATION 



A. General procedure and analysis 

Teleportation for continuous-variable systems [|[ is initiated when two parties, Alice and Bob, acquire two modes 
A and B, prepared in the joint quantum state pab, with joint Wigner function 



W PAB (a,p) = [^^tt[p AB D(a,ii)®D(b,>yj\D(j i ,a)D(>y ) 0) 



(3.1) 



Here Alice's mode is described by the annihilation operator a = (xa + ipA)/V% and Bob's by b = (xb + *Ps)/v / 2- 
We let a and (3 denote the respective complex-number variables for the modes. High-fidelity teleportation depends 
critically on the entanglement contained in pab- Thus the preparation of pab is generally done at some central point, 
after which the two modes are distributed to Alice and Bob. 

The state to be teleported is held by a third party called Victor. He brings up to Alice a mode V with annihilation 
operator v — (xy + ipv)/V^', the corresponding complex amplitude is denoted by v. Victor's mode is prepared 
in a state p, with Wigner function W p (u). For the present, we allow p to be pure or mixed; we specialize to pure 
input states when we introduce the teleportation fidelity below. The overall Wigner function for the three modes is 
W p (u)W PAS (a,0). 

Alice now measures the commuting variables xy + xa and py — pa, or more succinctly, she measures the complex 
quantity 

1 i 

w + a T = — (x v +x A ) + 7/| (PV ~Pa) ■ (3.2) 

The outcome, denoted by the complex number £ = (£i+i£2)/v2> occurs with probability p{£,)d 2 $,, where the probability 
density is given by 

p(0 = J ' d 2 v d 2 ad 2 f36(is + a* - £)W p (v)W PAB (a,p) = J ' d 2 vW p (v)W PA {€ - v*) (3.3) 

and pa — ^bPab is the marginal state of Alice's mode. The state of Bob's mode after the measurement, denoted by 
/?'(£), is conditioned on result £ and has Wigner function 

W p ,{(3\§ = -jrpjj d 2 vd 2 a8(v + a* -Z)W p {v)W PAB {a,fi) = J d 2 uW p (v)W PAB {C ~ v* , P) ■ (3.4) 

Alice now sends the measurement result £ to Bob, who displaces his mode by this amount; i.e., xb is displaced by 
£i and pb is displaced by £2, giving an output state Pout(£) with Wigner function 



W Po M£) = Wstf-ftt) = — / d*uW p (v)W PAB (C-S,l3-Z). (3.5) 



1 

pTO 

The average of the output state over the possible measurement outcomes 



Pout = J d^p(OPout(0 , (3.6) 

has Wigner function 

W Po J(3) = J ' d 2 £p(OW Pau MO = J d 2 V P(v)W p (f3- V ) = j ' d 2 vP{fS-v)W p {v) , (3.7) 

where 

P(u) = J d 2 aW PAB (a,v- a*) = j ' d 2 a d 2 (3 8(13 + a* - v)W PAB (a, (3) (3.8) 

is the probability density for obtaining result v in a measurement of the commuting observables xb +xa and ps ~Pa, 
i.e., in a measurement of 

1 i 

b + a' = -j=(x B + x A ) + ^{pb - Pa) ■ (3.9) 



5 



Since W p {(3—v) is the Wigner function for the displaced state D(v, v)pD" s (v, v), i.e., W p (/3—v) = W D r VtU \ pD i i VtV ){P), 
Eq. (|3.7[) implies that the average output state (|3.6p can be written as 

p out = I ' d 2 vP(v)D{b,v) P D\b,v) , (3.10) 



where in this equation we regard the initial state p as a state of Bob's mode. The average output state is an average 
over displaced input states, the average controlled by the distribution P(v). To achieve high fidelity in teleporting p, 
no matter how it is oriented in phase space, we need P[y ) to be a narrow distribution in all directions in phase space, 
highly peaked at v = 0. Throughout this paper, what we mean by high-fidelity teleportation is this ability faithfully 
to teleport p and any rotation of p. 

The symmetrically ordered characteristic function for p out is 

$p ou » = tr[p out D(6, M )] = 7rP0u)$», (3.11) 

where 

P(p) = J d ^P(y)D{v,p) (3.12) 

is the Fourier transform oiP{y). Equation (|3.11|) is simply the Fourier transform of the corresponding Wigner- function 
relation (|3.7| . We emphasize that P(p) is not a normalized probability distribution; rather its important properties 
are that ir\P(p)\ < 1 and irP(0) = 1. Teleportation with high fidelity requires that irP{p) be close to 1 over the entire 
region for which <& p (p) is nonnegligible. 

Up to this point, we have not needed to say whether the input state is pure or mixed, but we now want to measure 
the success of teleportation in terms of the overlap of the input state with the output state, called the fidelity. For this 
purpose, we need to assume that the input state p = is pure. We maintain this assumption until Sec. lVIl where 

we generalize our results to mixed states by using the entanglement fidelity in place of the fidelity. For outcome £, 
the fidelity of the output state with the input state is defined to be -F p (£) = (fp\Pout(0\''P) ■ Thus the average fidelity 
between input and output is 



F p = J d^p(0F p (0 = (VlTVtlV') • (3.13) 
We can manipulate the average fidelity into several forms, all of which play a role in the subsequent discussion: 

F p = / d 2 v P{v)\{^\D{b,v)\^)\ 2 = J d 2 vP{ V )\<5> p { V )\ 2 (3.14a) 
d 2 (id 2 vP{fi- v)W p {l3)Wp{v) (3.14b) 
ttI d 2 (3W Po J/3)W p (P) = irf d 2 l3d 2 vP{l3-v)W p {P)W p {v) (3.14c) 



7T 



The first line, Eq. (|3.14aj) . follows directly from inserting the average output state (|3.10p into the expression (|3.13p 
for the average fidelity. The second line, Eq. (|3.14b[) . comes from Fourier transforming the quantities in the integrand 
of the first line. The third line, Eq. (|3.14c[) . comes from rewriting (V'lPoutlV') as an overlap of the input and output 
Wigner functions, W Poxit (/3) and W P {P), and then using Eq. ([3~7| for W Pout (f3). Similarly, the fourth line, Eq. (|3.14d|) . 
comes from rewriting ("0|p out |'0) as an overlap of the input and output characteristic functions, <!>p out (/x) and $ p (p), 
and then using Eq. (|3. 1 1|) for Qp out (/i). The third and fourth lines are related to one another by a Fourier transform 
of the quantities in the integrand. 

These forms for the teleportation fidelity can be related to Zurek's work on sub-Planck structures in phase space • 
For a given input state p, the Wigner function W p (v) has two important length scales: (i) a small scale I over which 
the Wigner function varies substantially and (ii) a large scale L, which is the extent of the region over which the 
Wigner function is nonnegligible. In the characteristic function & p (p), these scales appear inversely: (i) a small scale 
7r/L over which, in some phase-space direction(s), the characteristic function plunges from 1 at p, = to close to zero 
and (ii) a large scale it /£, which is the extent of the region over which the characteristic function remains nonnegligible. 
As pointed out by Zurek [l6|, for pure states the large and small scales are reciprocally related, i.e., IL ~ 1. We make 
this reciprocal relationship explicit in Sec. HVl 
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In Eqs. (|3.14p . the first and second lines form a pair under Fourier transformation of the quantities in the integrand, 
and the third and fourth lines constitute another such Fourier pair. The first and second lines relate the average 
fidelity to the large-scale extent of the Wigner function W p (v), since good fidelity in Eq. (|3.14b|) requires irP to be 
close to 1 over the entire extent of the Wigner function. The third and fourth lines relate the average fidelity to 
the fine-scale structure of the Wigner function, since good fidelity in Eq. (|3. 14c|) requires P to be narrow relative to 
the fine-scale structure in the Wigner function. The strict connection between the two pairs and thus between the 
large-scale and small-scale properties of the Wigner function comes ultimately from the ability to write the average 
fidelity either as an overlap or by using the expression (|3.10l) for the average output state. 



B. Squeezed-state teleportation 

We now specialize the above analysis to the case of squeezed-state teleportation, where Alice and Bob choose a 
two- mode squeezed state for pab- We do a good job of teleporting when the distribution P(v) is narrow, i.e., when 
xa and Xb are tightly anti-correlated and pa and ps are tightly correlated. Introducing the modes c — (a + 6)/v2 = 
(xc + ipc)/V% and d = (a — b)/\/2 = (xd + ipD-)/\/2, corresponding to c-number variables 7 = (71 + i'j2)/V2 
and S = (Si + 162) /V%, respectively, we see from Eq. (|3.9j) that we want the variances of xc = (xa + xb)/V% and 
Pd = (pa — Pb)/ %/2 to be small. Thus a natural choice is to use squeezed vacuum states for modes c and d, with the 
variances of the quadrature components given by 

A 2 1 -2r A 2 1 2r 

Ax c = 2 e > A Pc = 2 e ' 

Ax D 2 = l -e 2r , A PD 2 = l -e~ 2r , (3.15) 

where r is the standard squeezing parameter. This state is a two-mode squeezed vacuum state for modes a and b with 
Wigner function 



W PCD (l,S) = ^expi-e^l-e-^i-e-^Sj-e 2 ^ 2 



= -2exp(-2|/3 + a*| 2 A-i|/?- a *| 2 /2) = W PAB (a,f3), (3.16) 

where we use the correspondence 7 = (a + /3)/\/2 and 5 — (a — /3)/v2 and introduce a new squeezing parameter 
t= 2e~ 2r . 

Since P(v) is the probability density to obtain result vina measurement of b + a^ — xc ~ipD, the Wigner function 
immediately implies that 

P(v) = A e - 2 M 2 /* : (3.i 7 ) 



with Fourier transform 



P(p) = I e -'H 2 / 2 . (3.18) 



Notice that the two- mode Wigner function (|3.16[) can be written as the product of these broad and narrow Gaussians, 
i.e., W PAB (a, 13) — 2tP((3 + a*)P((3 — a*). The marginal Wigner function for mode A, 

/2t f 2t \ 

d 2 f]W PAB (a,P) = 7T{1 + t2/A) ex P(^TT7V4 H2 J ~ 2<P ( 2a )' ( 3 - 19 ) 

is always broader than P(2a), but not by much for large squeezing. The final approximation holds in the limit of 
large squeezing. 

We can now specialize the important results of the preceding analysis to the case of squeezed-state teleportation. 
The average output state p.iop at Bob's end becomes 

Pout = \ f d ^e' 2 ^^D(b^) P D\b^) . (3.20) 
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With s — —t, the symmetrically ordered characteristic function for p out becomes the s-ordered characteristic function 
for p, 

$p ou » = vrP^pQi) = e-'W 2 / 2 ^) = (3.21) 
and hence the Wigner function for p ou t becomes the s-ordered quasidistribution for p, 

W~ Pa JP) = W^\P) = \ J ^e-^-^/^WM . (3.22) 
The various forms for the average fidelity in Eqs. (|3.14[) become 

F P (t) = \ f^e-^l^ p {u)\ 2 (3.23a) 

= J d 2 f3d 2 is e-^-^^Wpi^Wpiu) . (3.23b) 

= irj d 2 pwj, s Hp)W p (P) = -( d 2 [3d 2 ve - 2 ^^ t W p (f3)W p (v) (3.23c) 

= / V^'^^W = / ^e-^ 2 ' 2 \<S>M\ 2 (3.23d) 

The first and fourth lines (or the second and third) show that the teleportation fidelity obeys a scaling relation: 
tF p (t)/2 = F p (A/t). It is worth noting that these forms of the squeezed-state teleportation fidelity show that it is 
invariant under phase-space displacements and rotations of p. 

The relevant range of t is < t < 2. When t = 2, Alice's and Bob's modes are unentangled, each being in 
the vacuum state. In this case, the teleportation process reduces to a "classical" version of teleportation. Alice's 
joint measurement on her mode and Victor's mode is a heterodyne measurement, for which the probability density 
to obtain outcome £ is given by the Husimi function of p, i.e., p(£) = Q p (£,)- When Bob receives the outcome £, 
he displaces his vacuum state by £, yielding as output the coherent state p O ut{0 = The resulting fidelity, 

F p (£) = (£|p|£) = 7rQp(£)j is also given by the Husimi function, yielding an average fidelity 

F p {2) = nJd 2 ZQ 2 p (t) = J = / ^ e-M 2 \^> p ( P )\ 2 , (3.24) 

in agreement with Eq. (|3.23d[) . Values t > 2 are not unphysical, but they do correspond to antisqueezing of the 
quadratures xc and po and thus to fidelities worse than this classical process. 

C. High-fidelity limit and communication requirements 

In the limit of large squeezing, with r going to infinity and t going to zero, continuous-variable teleportation 
achieves high fidelity. The function P(v) becomes a very narrow Gaussian, approximating a delta function, while 
2tP(p.) becomes a very broad Gaussian. The average output state (|3.20p is obtained from the input state by Gaussian 
phase-space displacements of characteristic size \Jt/2 = e~ r '. The first form of the average fidelity, Eq. (|3.23a|) . shows 
that to get high fidelity between input and output, e~ r should be somewhat smaller than the size ir/L of fine-scale 
structure in the characteristic function of p; the fourth form, Eq. (|3.23d[) . assures us that the fidelity is near one when 
e r is somewhat larger than the extent tt/£ of the characteristic function. Thus the first and fourth forms express the 
reciprocal relationship between the fine-scale structure and large-scale extent of a pure-state Wigner function. 

These same conclusions can be read off the Wigner-function forms of the fidelity. The second form of the average 
fidelity, Eq. (|3.23b|) . tells us that the fidelity is close to 1 when e r is somewhat larger than the extent L of the Wigner 
function of p, whereas the third form, Eq. (|3.23c[) , assures us of high fidelity when e~ r is somewhat smaller than 
the scale I of fine-scale structure in the Wigner function. High fidelity is achieved when the available squeezing is 
sufficient to telcport faithfully the fine-scale structure in the Wigner function. 

In the high-fidelity limit, we can simplify the account of the teleportation process. To see what is going on, we take 
a closer look at the probability density p(£) for Alice to get outcome £ in her measurement and at the conditional 
Wigner function Wp out (/3|£) of the teleported state. The Wigner function (|3.19p of Alice's marginal state is a broad 
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Gaussian, and thus the probability density of Eq. I|3.3[) reduces to the same broad Gaussian, displaced to account 
for the location of W p {v), i.e., 

r _ Ae~ 2r 

p(0 ~ 2i \ d 2 vW p (y)P(2t-2v) ~ 2tP{2(_ - 2{v)) = — exp (-4e~ 2r |£ - («>| 2 ) ■ (3-25) 

In the second step, we take advantage of the fact that the broad Gaussian P is nearly constant over the extent of 
W p (v) to evaluate v in the broad Gaussian at a typical point within the extent, taken here to be the mean value 
v = (v) of Victor's mode. The Wigner function (|3.5p of the output state then becomes 

2/ r 

W Pm MO = d 2 V W p {v)P{p- V )P{P + v-2i) 

~ J d 2 vW p {v)P([3-v) 

= w p M = v/ ^^WH^-H 2 ) • ( 3 - 26 ) 

In the second step we set j3 = v = (v) in the broad Gaussian P and use the above approximation for p(£). The 
result is that in the limit of high-fidelity teleportation, the output state pout(C) is independent of the measurement 
result £ and thus is the same as /? out . In the limit t — > 0, P(f3 — v) becomes a S function, giving perfect fidelity, i.e., 

Pout(£) = Pout P- 

In the high-fidelity limit expressed by these equations, we can give a a very simple Heisenberg-picture account 
of teleportation. The mean value of the measurement result £ is given by (£) = (v) , with corresponding variances 
A£i 2 = A^ 2 = e 2r /4. Following Alice's measurement, we know that v + = £, implying that b = — + xq — ipD = 
v — £ + — «Pd- Displacing 6 by £ then gives b — v + xc — ipD- Thus, after the teleportation process is completed, 
Bob's mode is identical to Victor's initial mode, except for contamination by the fluctuations in xc and pd- Since 
the variances of xc and pr> are both equal to e~ 2r /2, the convolution (|3.26[) describes the output state corresponding 
to this overall transformation of Bob's mode. 

If Bob is to take full advantage of the squeezing resource, the value of £ transmitted to him must allow him to 
perform displacements with accuracy somewhat better than the uncertainties in xq and pu ■ Thus the required number 
of bits in the transmission of £i or £2 must be roughly log(A£i/Aa;c) = log(e 2r /\/2), making the total amount of 
information in transmitting both £1 and £2 approximately equal to 21og(e 2 ' /v / 2) — 4r/ln2 bits. The squeezing 
parameter must be large enough to teleport the smallest phase-space structures faithfully; i.e., e~ r must be somewhat 
smaller than the fine-scale structure in the Wigner function of p, of size t ~ 1/L. Thus Alice must transmit roughly 
2 log L 2 bits in order to achieve high-fidelity teleportation. 

Since L 2 is approximately the number of phase-space-local states needed to represent p, the number of classical bits 
transmitted has the standard form of 2 bits for each qubit of quantum information. Indeed, this argument gives us 
an independent way of interpreting the required 2 bits of classical information per qubit. The teleportation process 
must be able to distinguish (Ljt) 2 phase-space regions to transmit all of the sub-Planck structure, and this means 
transmitting log(L/^) 2 ~ 21ogL 2 bits of classical information. 



IV. TELEPORTATION FIDELITY AND SUB-PLANCK STRUCTURE 

The discussion of high-fidelity teleportation in Sec. IIII CI draws attention to how the expressions (|3.23p connect 
teleportation fidelity to the fine-scale structure and the large-scale extent of the Wigner function of the state being 
teleported. In this section we make these connections explicit. 

We begin by noting that the first two derivatives of the average fidelity, as expressed in Eq. (|3.23d[) , 

^ = -\J d ^\tfe^ 2 ^ P {p)\ 2 < 0, (4.1) 
^ = i/^H 4 e- tM2/2 |$pM| 2 > 0, (4.2) 

imply that F p (i) is a strictly decreasing, strictly concave function oft. We can construct a characteristic scale for t at 
which teleportation becomes ineffective by approximating the average fidelity for small t as F p (t) ~ 1 + t(d F p /dt)t=o 
and asking when this approximation goes to zero. The result is a critical value of t given by t c — | (d F p j dt)t=o | 1 ■ For 
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t > t c , teleportation becomes ineffective because the avaiiabie squeezing is unable to resolve the fine-scale phase-space 
structure of p. Since t — 2e~ 2r , we convert this critical value into a phase-space length that characterizes the size of 
the fine-scale structure by defining 



£ c = y/tc/2 



1 



(4.3) 



>2\(dF p /dt)t=o\ 

To evaluate £ c , we manipulate the derivative (|4.ip . evaluated at t = 0, through the following sequence of steps: 



dF n 



dt 



t=o 



\ J ^Vl* P &«)| s 



^d 2 08^-fx)fx*P^)^) 



~ d 



Noticing that 



we can write 



dh 



dW n 



dW p {v) 



dv 



1 f dW p ,dW P 



s/2 \ dv x 



dv 2 



dW n 



(dW p 



dW n 



= l\VW p ( Vll v 2 )\ 2 



Putting this together, we find that 



j J du x dv 2 | VW p (vi , v 2 ) 



-1/2 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



This measure is motivated here by the decrease of teleportation fidelity, but it is a reasonable a priori measure of 
the linear size of fine-scale structure in the Wigner function. As noted in the Introduction, for a general normalized 
distribution or quasidistribution, £ c quantifies the size of the fine-scale structure multiplied by the square root of the 
area of support of the distribution. Since all pure-state Wigner functions have the same area of support, corresponding 
to one Planck area, £ c can be interpreted as measuring of the size of the fine-scale structure. 

We can also evaluate £ c using the upper Fourier pair in Eqs. (|3.23[) . and this will relate our measure of fine-scale 
structure to the large-scale extent of the Wigner function. Differentiating and manipulating Eq. (13.23bj) . we have 



dF 



t=o 



d 2 f3d 2 is\(3-is\ 2 W p (l3)W p (v) 



■t / d 2 pd 2 v (|/3| 2 + M 2 - Pv* - P*y)W p ((i)W p (v) 



l((vh) + (vvi)) + \(v)\* 



which gives 



where we define 



2 

-- (Ax v 2 + A Pv 2 ) 



£ c = (Ax v 2 + Ap< 



2\-!/2 



2 

T,. 



L c = 2\/Ax v 2 + A Pv 2 



(4.8) 



(4.9) 



(4.10) 
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as a sensible measure of the linear extent of the region over which the Wigner function is nonnegligible. Through the 
use of the subscript c, we emphasize the explicit definitions of £ c and L C1 as opposed to the somewhat loose use of I 
and L up until now. The relation 

i c L c = 2 , (4.11) 

arising from our considerations of teleportation fidelity, is a rigorous expression of the reciprocal size of the fine-scale 
and large-scale structures of a pure-state Wigner function. 

It is easy to see that Axy 2 + Apy 2 = (Axy — Apv) 2 + 2A.xyIS.py < 2AxyApy < 1) where the last step is the 
Heisenberg uncertainty principle. Equality holds in both inequalities if and only if p is a coherent state. Thus coherent 
states have the smallest initial rate of decrease of average fidelity as t increases from zero. Equivalently, they have the 
smallest extent, L c — 2, and no fine-scale structure, i.e., £ c = 1. That coherent states have the smallest initial rate of 
decrease of fidelity is a reflection of the fact that they have the highest teleportation fidelity for all values of t, a fact 
we demonstrate in the next subsection. 



V. EXAMPLE STATES 



We now examine various simple examples to develop our intuition for the relationship between teleportation fidelity 
and sub-Planck structure. In the following we investigate coherent, squeezed, and number states, Zurek's compass 
state, a class of random states, and a time-evolved state of a chaotic system. 



A. Coherent states 



All coherent states give the same average fidelity as the vacuum state, for which 

$|o>0l0 = (0|£>(M)|0) = e-^ 2 / 2 . (5.1) 
Therefore, the teleportation fidelity for any coherent state \v) is 

r d M r -(i+t/2)\n\ 2 _ 

- J n 

and 4 = 1 (L c = 2). 

Coherent states achieve the maximum teleportation fidelity for all values of t, a fact we pause to demonstrate. To 
do so, we notice that the average fidelity in the form (|3.23b[) can be thought of as the average value of e~*'^~^ I 2 
with respect to a product state p ® p of two modes, b and v, the joint Wigner function for the two modes being 
Wbv(P, v) — Wp{j3)W p {v). What we actually find is the maximum value of this average value for all two-mode states 
Pbv i n ot just tensor products of copies. 

Using the fact that the Wigner function returns expectation values of symmetrically ordered operator products, we 
can write the relevant average as the expectation value of an operator A t , 

d 2 f3d 2 ve- l ^ 2/2 W PBV {f3,v) = tr( PBV A t ), (5.3) 



where 

l + t/2\l + t/2j ' V ; 

with / = (b-v)/V2. The maximum of the expectation value is given by the largest eigenvalue of At ■ The eigenstates 
of A t are the number eigenstates of the mode with annihilation operator / in tensor product with any state of the 
mode with annihilation operator g = (b + v)/y2. Since the factor in large parentheses in the expression (|5.4p has 
magnitude < 1, the largest eigenvalue is (1 + t/2)^ 1 , which occurs (uniquely when t > 0) for the vacuum state of 
mode /. This establishes that the maximum teleportation fidelity is the coherent-state fidelity (|5.2[) . 

By returning to the case of interest, i.e., psv being a tensor product of two copies of a pure state, it is easy to show 
that for t > 0, the fidelity bound is saturated only by coherent states. 
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FIG. 1: Zurek's compass state with a = 5/\/2 as depicted by (a) |$ P | 2 , (b) (ty/2)W p , and (c) irQ p . The Wigner function displays 
fine-scale structure on the scale £ c = l/\/2a = 0.20, and the size of its large-scale extent is given roughly by L c — 2/£ c — 10. 
These scales are reversed in the characteristic function, which has an extent characterized by n/£ c — nL c /2 — 16 and fine- 
scale structure on a scale i:/L c = n£ c /2 — 0.31. In (a) and (b), the scale bars indicate ttL c /2 and L c , respectively, and the 
insets, which are blow-ups of the regions bounded by dashed lines, show the fine-scale structure in more detail, with scale bars 
indicating n£ c /2 and £ c , respectively. 



B. Squeezed states 

All squeezed states with the same squeezing parameter u have the same average fidelity; so we need only calculate 
the fidelity for a squeezed vacuum state 



|0,«) = e"( &2 - bt2 )/ 2 |0) . 

Using 

e -u(b 2 -b^)/2 D ^ b ^yu(b 2 -b^)/2 = D(b,ticoshu + ii*sinhu) = D(b, (e>i + ie-> 2 )/V2) 
we find the characteristic function, 

$|0,»)(M) = e -( e2 >?+^ 2 >2)/ 4 . 

The teleportation fidelity for any squeezed state with squeezing parameter u is then 

72„ ....... 1 



F\is,u)(t) 



d fj e -t|p| 2 /2 e -(e 2 >?+e- 2 "^)/2 



y/l + tcosh2u + t 2 /4 



(5.5) 



(5.6) 



(5.7) 



(5.8) 



giving £ c = l/\/cosh2u (L c — 2\/cosh2u). For a highly squeezed state (u^> 1), we have £ c ~ y/2e u and L c ~ y^2e 



C. Number states 



The number-state matrix elements of the displacement operator are given by [25j 

lm\D(b u)\n) fV^-e- M2/2 ^- n Lt- n) m, i£m>n, 



(5.9) 



where L™{x) are generalized Laguerre polynomials [26| . The characteristic function of a number state is a diagonal 
matrix element: 



* N 0*) = (n\D(b,fx)\n) = e-^ 2 i 2 L n (\tf) 



(5.10) 
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FIG. 2: (a) Average fidelity (|5.16[1 of the compass states for a = 0, 0.1, 0.2, . . . , 5 (lighter tones from top to bottom) and of the 
compass state in Fig. Q] (full dark). The upper and lower bounds from Eq. (|5,16[1 are also shown (dashed), as is the level of 
fidelity at the critical squeezing parameter t c — 2£^ (dash-dotted). When a — the compass state is a single coherent state 
centered at the origin with fidelity (upper dashed curve) given by Eq. (|5.2|l . The a — + oo bound 1/4(1 + t/2) from Eq. (15.161) 
is shown as the lower dashed curve, most of which is obscured by the full dark line for the state of Fig. [1] The dash-dotted 
curve shows that to achieve a fidelity of approximately 1/2 or better, we need to teleport with a squeezing parameter t < t c . 
(b) Average fidelity of the compass states as a function of a for t — 0.2,0.4, ... ,6 (lighter tones from top to bottom). The 
fidelity declines sharply when a = y/ir/2 (dashed), at which point the four coherent states are separated by a distance specified 
by a von Neumann lattice. This is the separation at which the interference fringes and checkerboard pattern of Fig. [T] appear. 



The number-state fidelity is thus 



F„ (/-. = / ^e-^ 2 ^[L n (\^)] 2 (5.11) 
dx e-^ +t ' 2 >[L n {x)] 2 (5.12) 
^f(l + i/2)- 2 »- 1 J F( - n, -n; -2n; 1 - t 2 /A) (5.13) 



(i-*/2)» (i+e/A\ 

(l+t/2) n + 1 n \l-t 2 /4J 



2.„ , 

(5.14) 



where F(a, b; c: x) is a hypergeometric function, P n (x) is the Legendre polynomial, and we have used formulae found 
in Refs. [26|, [231 . Using Eq. (|4.3p . we find that l c — 1/ ' \j2n + 1. More directly, we can use the number-state variances 
(Ax) 2 = (Ap) 2 = n + \ to find 2/4 =L C = 2 v / 2rT+T. 



D. Zurek's compass state 

Zurek introduced the "compass state" in his original article on sub-Planck structures It is a superposition of 
four coherent states at positions (x,p) — (\/2a,0), (—^2^,0), (0, V^a), and (0, — \/2a): 

m = k) + \-a) + H + Ha) , 5 15) 

2e- a2 /2^/2(cosha 2 + cosa 2 ) ' 
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FIG. 3: A random state with N = 100 as depicted by (a) |3v| 2 , (b) (tv/2)W p , and (c) ttQ p . The Wigner function displays fine- 
scale structure on the scale l c = 1/vN = 0.10, and the size of its large-scale extent is given approximately by L c = 2/£ c = 20. 
These scales are reversed in the characteristic function, which has an extent characterized by ir/£ c = nL c /2 = 31 and fine-scale 
structure on a scale ir/L c — ir£ c /2 = 0.16. The scale bars and insets in (a) and (b) are as in Fig. [1] 



A straightforward, but laborious calculation yields the average teleportation fidelity: 



*]*>(*) 



4(1 + t/2) 



2 £ ^ 2 t o i _ ( n 2 i p \ ( o 

cosh a + cos a +2 cosha + cos a cos a + cosh a 

2+t 2+t J \ 2+t J\ 2+t 



(cosha 2 + cos a 2 )' 



l + t/2 ' 
1 

4(1 -M/2) 



fora->0, 
, for a — > oo . 



(5.16) 



Differentiating this expression gives 



L = 1 + 2a 



, sinh a 2 — sin a 2 
cosh a 2 + cos a 2 



-1/2 



y/2a 



for a 



(5.17) 



More directly, one can evaluate (Ax) 2 + (Ap) 2 and use Eq. (|4.8|) to find ^ c = 2/L c . 

The Wigner function of a compass state with a = 5/ \/2, which is well into the regime of large a, is shown in Fig. [TJb) . 
Interference fringes form between adjacent coherent states, combining at the origin to create a checkerboard of fine 
structure at a scale on the order of £ c = \j\f2a = 1/5. Complementary behavior is displayed in the Fourier transform, 
the characteristic function, whose absolute square is plotted in Fig. QJa). The Husimi function [Fig. [ljc)] can be 
viewed as a Gaussian smoothing of the Wigner function. The average teleportation fidelity (|5.16[) for a compass state 
is plotted in Fig. [Ha) for various values of a. The caption explains in detail the various features plotted in Fig. [51 



E. Random states 



We now consider states of the form 

N-l 

|V) = ^cn|n) (5.18) 

where the states \n) are the number eigenstates and the coefficients c n form a a random complex unit vector in N 
dimensions under the uniform measure, i.e., a random point on the unit sphere in 2N dimensions. These states are 
conjectured to have the same statistical properties as an eigenstate (or long-time evolved state) of a chaotic system 
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FIG. 4: The average fidelity (|5.24|l for random states with N = 1, 2, ... , 100 (lighter tones from top to bottom). The case N — 1, 
which reduces to the coherent-state fidelity (|5.2|l . is highlighted (dashed), and the level of fidelity at the critical squeezing value 
t c is plotted (dash-dotted), showing again that to achieve a fidelity of approximately 1/2 or better requires t < t c . The fidelities 
of the particular random state in Fig. [3] (dark full) and the chaotic state in Fig. [S] (dotted) are also drawn; they are essentially 
indistinguishable . 



with a region of ergodicity within a circle of radius r = y/2N in phase space [HI]. Figure [3] plots the absolute square of 
the characteristic function, the Wigner function, and the Husimi function for an example random state with N — 100. 

To calculate the average fidelity of a random state, we use the fidelity in the form (|3.23d[) , noting that the absolute 
square of the characteristic function for any state of the form fj5 . 18[) is given by 



JV-l 

|$ P ( M )| 2 = MD(b,^)\ 2 = ]T c^CjcUrnlDib, (1)1x1) (k\D^b^)\j) , 

m,n,k,j—0 

The necessary averages over random states as described above are given by 

-ppr * sjc-i ^mn^jk ^mj^nk 

HCnC m C 3 C k \ = N{N+l) ■ 

These averages imply that E[c„c*„] = S nm /N. As a result, we have 

1 N ~ 1 



m,n—0 



Using Eq. (|5.9[) and Ref. [27| . one can show that 
-e - tM2/2 (m\D(b, [i)\ri)(n\Di(b, fi)\m) 



(to + n)\ 



TV 

;2, 



< L-H e - t ^ 2/2 {m\D{b, ii)\m)(n\DUb,iJL)\n) 

TV 



m\n\ (1 + t/2) m+n + 
(m + n)\(t/2) m+n 



which gives 



N-l 



E 



N(N + 1) m! n\ (1 + i/2) ro +™ +1 



(m + n) ! 



m!n! (1 + 1/2) 



— to, — n; —to — n; 1 — ^ 2 /4) 



Y-F( — to, — n; — m — n; 1 — i 2 /4) 
— ™i — n ; — m — t&; 1 — 4/i ) 



+ (t/2) m+n F( - to, -n; -m - n; 1 - 4/i 2 ) 
This average fidelity is plotted in Fig. @]for iV = 1,2, ... , 100 (lighter tones). 
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FIG. 5: A chaotic state produced by the driven double-well Hamiltonian (|5.28[) as depicted by (a) |"Iv| 2 , (b) (ir/2)W p , and 

(c) TYQp. 



To calculate the slope of the fidelity at t = 0, it is easiest to work from Eq. (|4.8j) to obtain 

N—l N-2 

2 



N-l N-2 

c* m+1 . (5.25) 



t=o 



Averaging over random states leads to 
1. 



-E[(Ax) 2 + (Ap) 2 ] = E 







dt 


t=0- 



m,n— 



dE[F m ] 



dt 



We now define our measures of small- and large-scale structure to be 



2/4 = L c = 2VE[(Az) 2 + (Ap) 2 ] = 2 



A^ 2 + 1 
N + 1 



A" 2 + 1 
2(JV+ 1) 



2VA^ for A^ -> oo , 



(5.26) 



(5.27) 



ignoring that the averaging over random states doesn't commute with taking square roots and reciprocals, although 
it will do so very closely when Af is large. 

F. Chaotic state 

We now consider briefly sub-Planck structure in states of a chaotic Hamiltonian. Consider a long-time evolved state 
of the driven double-well potential, corresponding to Hamiltonian 



H = 5p 2 - 8.t 2 + 0.05x 4 + 65x cos(27rf) 



(5.28) 



For the parameters chosen in this Hamiltonian, this system is chaotic in the classical limit. By choosing an initial 
coherent state centered at (x,p) = (—8,4) and then evolving the system from t = to t = 5, we obtain the chaotic 
state pictured in Fig. [5l The Wigner function displays an abundant fine-scale structure that qualitatively resembles 
that of the random state depicted in Fig.[3l It is no surprise that the fidelity curve also follows that for random states. 
This is plotted in Fig. [4] as the dotted curve, but is obscured behind the fidelity for the random state (full dark). 

VI. MIXED-STATE TELEPORTATION AND ENTANGLEMENT FIDELITY 



In the preceding sections the teleported state was assumed to be pure. Since teleportation is a linear operation, 
the procedure outlined in Sec. IIIII works equally well for mixed states. The overlap between input and output states, 
however, is no longer an appropriate measure of teleportation fidelity. A suitable measure for assessing the fidelity 
with which a mixed state p is teleported is the entanglement fidelity [23| |. which is the fidelity for teleporting Victor's 
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half of a purification of p, thus transferring the entanglement to Bob (entanglement swapping). It is quite easy to 
see how to generalize all of our results to entanglement fidelity. Victor's mode is now entangled with another mode, 
labeled by U, which has annihilation operator u and corresponding complex variable p. The joint state of U and V is 
a pure state puv = \i>uv)('4>uv \, which purifies Victor's state, i.e., tru(puv) = P- The results of Sec. [Ill] [Eqs. (|3.3D . 
(|3.5|) . (|3.7|) . (|3.10p and (|3.11|) ] generalize immediately to the following: 



W Pout (»,f3\® 



d 2 vW puv {tx,v)W pAB {e , 



d 2 pd 2 vW Puv (p,v)W PA {C -v*) 

1 

W Pmit (n,(3) = I ' d\P{v)W Pvv {p^-v) 

Pout = 



*p ou >, a) 



d 2 vP{v) [I ® D(b, v)] puv [I <S> D'(b, u)] , 
■nP(a)$ Puv {v,a) , 



The entanglement fidelity averaged over outcomes £ is 



= WwlTwIVw) = ^ / d 2 pd 2 pw- Pmt (p,p)w PU A^P) 



(6.1) 

(6.2) 

(6.3) 

(6.4) 
(6.5) 

(6.6) 



notice that when p is pure, so that \ipuv) is a product state, the entanglement fidelity reduces to the ordinary 
fidelity (|3.13[) . The average entanglement fidelity can be put in the four forms analogous to Eqs. (|3.14p : 



d 2 vP{v)\$ p {v)\ 2 
= nj d 2 (3d 2 vP{f3-v)W p {(3)Wp{v) 
= it J d 2 fid 2 vP{fi~v) irj d 2 pW puv (p 1 /3)W puv (p 1 iy) 



d 2 aP(a) 



d 2 v , 

— \$ puv (v, a)Y 



(6.7a) 
(6.7b) 
(6.7c) 
(6.7d) 



The first line, Eq. (|6.7a[) . comes from inserting p out into the first form of the fidelity in Eq. (|6.6|) . Fourier transforming 
the terms in the integrand yields the second line, Eq. (|6.7b[) . This Fourier pair is identical to the upper two lines in the 
pure-state teleportation fidelity (|3.14p ; they relate the entanglement fidelity to the large-scale extent of the Wigner 
function of p or, equivalently, to the fine-scale structure of the characteristic function of p. The third and fourth lines 
come from writing the entanglement fidelity as an overlap of joint Wigner functions or joint characteristic functions. 
These two lines also constitute a Fourier pair, which relates the average entanglement fidelity to the fine-scale structure 
in the joint Wigner function W Puv (p, (3) or, equivalently, to the large-scale extent of the joint characteristic function 
$ Puv (y, a). This Fourier pair is different from the corresponding third and fourth lines for the pure-state fidelity (|3.14p 
precisely because these properties appear in the joint functions, not in the Wigner and characteristic functions for p 
alone. 

We can, however, convert the third and fourth lines to forms that involve only system operators. For this purpose, 
we focus first on the characteristic-function integral 



d 2 v 



d 2 v 



\$> puv {-v\a)t 



(6.8) 



The reason for the term on the right becomes clear as we manipulate this integral below. Since the entanglement 
fidelity is independent of which purification is used, we can use the purification 



\i>uv) = ^2 \n)u ® ^fp\n) v 



(6.9) 



n=0 



where \n)u and \n)v denote number states for U and V, respectively. 
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The characteristic function now becomes 

$ pc/v(-^*:«) = {^uv\D(u,-v*) <g> D(v,a) \ipuv) 

= ^2u(n\D(u, -v*)\m)uv{n\y/pD(v, a)^/p \m) v 

m,n 

= X! v(m\D(v, v)\n) v {n\^fpD{v, a)<Jp \m) v 

771, n 

= tT[D(v,v)y/pD(v,a)y/p] , 

where in the third step, we use 

u{n\D{u,-u*)\m)u = u {n\D(u 1 -v)\m)* u = u(m\D(u,i/)\n)u = v{m\D(v,i/)\n) v 
Inserting Eq. (|6.10p into the integral (|6.8[) gives us a form that involves only system operators: 



(6.10) 
(6.11) 



/ 



d 2 v _ _ 

\4> Puv (v,a)\ 2 = I — \tx[D{v,u)^pD{v,a)yfp) 

J 7T I 

f d 2 v 

= / ^v(m\^D(v 1 a) y /pD(v,iy)\m)v(n\D^(v,i')y/pD\v,a)y/p\n)v 

= ^2v(n\^/pD(v,a)^/p^/pD^(v,a)^/p\n) v 



d 2 v 



= tr[pD^(v, a)pD(v, a)] 
The simplification in the third step follows from 

d 2 v 



D(b, v)\m) v (n\D\b, v) 



(6.12) 

(6.13) 
(6.14) 

(6.15) 



which is a consequence of Schur's Lemma for the (Weyl-Heisenberg) group of displacement operators, but which can 
also be derived directly, for example, from Eq. (|5.9p . 

The joint Wigner-function integral in Eq. (I6.7cp can now be obtained by a Fourier transform of all variables in the 
characteristic- function integral (|6.14|) : 



tt/ d 2 (xW puv (n*,f3)W puv (n*,v) 



ti[pD(v,/3)pD(v,is)] = W p ((3,v) 



(6.16) 



Here D(v,v) is the Fourier transform of the displacement operator, as in Eq. (|2.6[) . and we define the two- variable 
function W p (/3, v). For a pure state p, W p (/3, v) = W p ((3)Wp(v) : and for mixed states, W p (/3, v) is what replaces the 
product of Wigner functions in the third form (|3.14c[) of the teleportation fidelity. The Fourier transform of W p ((3, v) 
is 



a) 



d 2 (3 d 2 vW p {J3,i>)D* {P, p)D(v,a) = tr[p£) t (v, fi)pD(v,a)] , 



(6.17) 



For a pure state p, we have & p (p, a) = &*(p)& p (a). 

We can now summarize our results by rewriting the third and fourth lines in the entanglement fidelity [Eqs. (|6.7c|) 
and Eqs. (|6.7d)) ]. the lines that tell us about fine-scale phase-space structure, in terms of the new bold-face functions: 



= n d 2 (3d z vP{(i-v)W p {(3,v) 



d 2 pP(p)$ p (p,p) 



(6.18a) 
(6.18b) 



These forms can be specialized to the case of squeezed-state teleportation by inserting the expressions for P and P 
from Eqs. (f3TTT|) and |3~T8)) . 

When we repeat the steps of Sec. lIVI to find the first derivative of the entanglement fidelity, the sequence of Eq. (|4.8j) 
is unchanged, because it uses the upper Fourier pair in Eq. (|6.7p . but the sequence of Eq. (|4.4[) . which uses the lower 
Fourier pair, must be modified to use the bold-face functions. The resulting expressions for the derivative, 



dF ent 



df 



(Ax v 2 + A Pv 2 ) 



t=Q 



d2y d 2 W p (v,a) 



dv* da 



(6.19) 
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give rise to the mixed-state generalizations of our measures for small- and large-scale phase space structure: 

V2 . , , . „,„, / \ ,1/2 



dF cnt 

-2 9 



dt 



1 



Lc 
2 



Ax v 2 + A Pv 2 



2v 8> WpM 



dv* da 



It is instructive to see how this works out for a thermal state, 



d 2 (3 



(6.20) 



(6.21) 



where A is the dimensionless inverse temperature and n = (e x — l)" 1 is the mean number of quanta. The second 
form is the standard P-function representation of a thermal state [251 ]. Inserting the P- function representation into 
the expression 1T[) gives 



1 



2n+ 1 



exp 



2n 2 + 2n+l,. ., . 2 , n(n + 1) 
- 2(2n + l) (W a + H a ) + -k+l 



(/ia* + 



and thus 



*(m,m) 



1 



-|M| 2 /(2n+l) 



2n+ 1 

The resulting average entanglement fidelity from Eq. ()6. 18b[) is 

1 



rp cnt 

p 1 + (2n + l)t/2 ' 



This gives 



(6.22) 



(6.23) 



(6.24) 



(6.25) 



which is consistent with the variances of x and p, i.e., (As) 2 + (Ap) 2 = 2n + 1. 



VII. CONCLUSION 



In this paper we examined the relationships among the output fidelity of continuous- variable teleportation protocols, 
sub-Planck structure in the Wigner function of the teleported state, and the large-scale extent of that Wigner function. 
For pure states, these relationships are made mathematically precise in Eqs. (|4.4[) and (|4.8|) . which lead us to define 
measures of small- and large-scale structure for the Wigner function in Eqs. (|4.7|) and (|4.10p . Consideration of 
several example states in Sec. [V] illuminates these relationships and builds confidence that the measures of small- and 
large-scale structure we define are quite reliable measures of phase-space properties of the Wigner function. 

The generalization of these results to mixed states in Sec. IVIII leads to a pair of new functions, generalizations of 
the Wigner and characteristic functions, which capture the fine-scale phase-space structure in any purification of the 
mixed state. These new functions might prove useful in other studies of phase-space properties of mixed states. 
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